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the Sun (obtained using long- baseline radio interferometry) , and ranging to Mars using the Viking 
lander, constrain the numerical values of the parameters of the specific models. 
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I. INTRODUCTION 

The idea, proposed in [1], that our four-dimensional 
Universe might be a three-brane, embedded in a five- 
dimensional space-time (the bulk), has attracted con- 
siderable interest in the past few years. According to 
the brane world scenario, the physical fields (electromag- 
netic, Yang-Mills, etc) in our four-dimensional Universe 
are confined to the three-brane. These fields are as- 
sumed to arise as fluctuations of branes in string the- 
ories. Only gravity can freely propagate in both the 
brane and bulk space-times, with the gravitational self- 
couplings not significantly modified. This model origi- 
nated from the study of a single 3-brane embedded in 
five dimensions, with the 5D metric given by ds^ = 
e~^^y^ri^^dx^dx'^ + dy"^, which due to the appearance of 
the warp factor /(y), could produce a large hierarchy be- 
tween the scale of particle physics and gravity. Even if 
the fifth dimension is uncompactified, standard 4D grav- 
ity is reproduced on the brane. Hence this model al- 
lows the presence of large, or even infinite non-compact 
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extra dimensions. Our brane is identified to a domain 
wall in a 5-dimensional anti-de Sitter space-time. Due to 
the correction terms coming from the extra dimensions, 
significant deviations from the Einstein theory occur in 
brane world models at very high energies [2, 3]. Gravity 
is largely modified at the electro-weak scale 1 TeV. The 
cosmological and astrophysical implications of the brane 
world theories have been extensively investigated in the 
literature [4]. Gravitational collapse can also produce 
high energies, with the five dimensional effects playing 
an important role in the formation of black holes [5] . 

In standard general relativity the unique exterior 
space-time of a spherically symmetric object is described 
by the Schwarzschild metric. In the five dimensional 
brane world models, the high energy corrections to the 
energy density, together with the Weyl stresses from bulk 
gravitons, imply that on the brane the exterior metric of 
a static star is no longer the Schwarzschild metric [6]. 
The presence of the Weyl stresses also means that the 
matching conditions do not have a unique solution on 
the brane; the knowledge of the five-dimensional Weyl 
tensor is needed as a minimum condition for uniqueness. 
Static, spherically symmetric exterior vacuum solutions 
of the brane world models were first proposed in [6] and 
in [7]. The first of these solutions, obtained in [6], has the 
mathematical form of the Reissner- Nordstrom solution of 
standard general relativity, in which a tidal Weyl param- 
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etcr plays the role of the electric charge of the general 
relativistic solution. The solution was obtained by im- 
posing the null energy condition on the 3-brane for a bulk 
having non zero Weyl curvature, and it can be matched 
to the interior solution corresponding to a constant den- 
sity brane world star. A second exterior solution, which 
also matches a constant density interior, was derived in 
[7]. 

Several classes of spherically symmetric solutions of 
the static gravitational field equations in the vacuum on 
the brane have been obtained in [8-11]. As a possible 
physical application of these solutions the behavior of 
the angular velocity vtg of the test particles in stable 
circular orbits has been considered [9-11]. The general 
form of the solution, together with two constants of in- 
tegration, uniquely determines the rotational velocity of 
the particle. In the limit of large radial distances, and 
for a particular set of values of the integration constants 
the angular velocity tends to a constant value. This be- 
havior is typical for massive particles (hydrogen clouds) 
outside galaxies, and is usually explained by postulating 
the existence of dark matter. The exact galactic metric, 
the dark radiation, the dark pressure and the lensing in 
the flat rotation curves region in the brane world sce- 
nario has been obtained in [11]. It is also interesting to 
note that the flat rotation curves can also be explained 
in f{R) modified theories of gravity, without the need of 
dark matter [12]. 

Furthermore, two families of analytic solutions of the 
spherically symmetric vacuum brane world model equa- 
tions (with gtt —1/grr), parameterized by the ADM 
mass and a PPN parameter /3 have been obtained in 
[13]. Non-singular black- hole solutions in the brane world 
model have been considered in [14], by relaxing the con- 
dition of the zero scalar curvature but retaining the null 
energy condition. The four-dimensional Gauss and Co- 
dazzi equations for an arbitrary static spherically sym- 
metric star in a Randall-Sundrum type II brane world 
have been completely solved on the brane in [15]. The 
on-brane boundary can be used to determine the full 5- 
dimensional space-time geometry. The procedure can be 
generalized to solid objects such as planets. A method 
to extend into the bulk asymptotically flat static spher- 
ically symmetric brane world metrics has been proposed 
in [16]. The exact integration of the field equations along 
the fifth coordinate was done by using the multipole (1/r) 
expansion. The results show that the shape of the hori- 
zon of the brane black hole solutions is very likely a flat 
"pancake" for astrophysical sources. The general solu- 
tion to the trace of the 4-dimensional Einstein equations 
for static, spherically symmetric configurations has been 
used as a basis for finding a general class of black hole 
metrics, containing one arbitrary function gu = A{r), 
which vanishes at some r = > (the horizon radius) 
in [17]. Under certain reasonable restrictions, black hole 
metrics are found, with or without matter. Depending 
on the boundary conditions the metrics can be asymptot- 
ically flat, or have any other prescribed asymptotic. For 



a review of the black hole properties and of the lensing 
in the brane world models see [18]. 

Now, to be viable models, the proposed models need 
to pass the astrophysical and cosmological observational 
tests. There are several possibilities of observation- 
ally testing the brane world models at an astrophys- 
ical/cosmological scale, such as using the time delay 
of gamma ray bursts [19] or by using the luminosity 
distance-redshift relation for supernovae at higher red- 
shifts [20]. 

In addition to these possibilities, we also mention work 
on the gravitational lensing in brane world models [21], 
on the role of the brane charge in orbital models of high- 
frequency quasiperiodic oscillations observed in neutron 
star binary systems [22], on the complete set of analytical 
solutions of the geodesic equation of massive test particles 
in higher dimensional spacetimes which can be applied 
to brane world models [23], and on brane world correc- 
tions to the charged rotating black holes, to the pertur- 
bations in the electromagnetic potential and test particle 
motion around brane world black holes [24], The clas- 
sical tests of general relativity, namely, light deflection, 
time delay and perihelion shift, have been analyzed, for 
gravitational theories with large non-compactified extra- 
dimensions, in the framework of the five-dimensional ex- 
tension of the Kaluza- Klein theory, using an analogue of 
the four-dimensional Schwarzschild metric in [25]. So- 
lar system data also imposes some strong constraints 
on Kaluza-Klein type theories. The existence of extra- 
dimensions and of the brane world models can also be 
tested via the gravitational radiation coming from pri- 
mordial black holes, with masses of the order of the lunar 
mass, M ~ 10~'^ Mq, which might have been produced 
when the temperature of the universe was around 1 TeV. 
If a significant fraction of the dark halo of our galaxy 
consists of these lunar mass black holes, a huge number 
of black hole binaries could exist. The detection of the 
gravitational waves from these binaries could confirm the 
existence of extra-dimensions [26] . 

It is the purpose of the present paper to consider the 
classical tests (perihelion precession, light bending and 
radar echo delay) of general relativity for static gravita- 
tional fields in the framework of brane world gravity. To 
do this we shall adopt for the geometry outside a com- 
pact, stellar type object (the Sun), specific static and 
spherically symmetric vacuum solutions in the context of 
brane worlds. As a first step in our study, we consider the 
classical tests of general relativity in arbitrary spherically 
symmetric spacetimes, and develop a general formalism 
that can be used for any given metric. In particular, 
we consider the motion of a particle (planet), and ana- 
lyze the perihelion precession, and in addition to this, by 
considering the motion of a photon, we study the bending 
of light by massive astrophysical objects and the radar 
echo delay, respectively. Existing data on light-bending 
around the Sun, using long-baseline radio interferome- 
try, ranging to Mars using the Viking lander, and the 
perihelion precession of Mercury, can all give significant 



3 



and detectable Solar System constraints, associated with 
the brane world vacuum solutions. More precisely, the 
study of the classical general relativistic tests, constrain 
the parameters of the various solutions analyzed. 

This paper is organized in the following manner. In 
Section II, we outline, for self-completeness and self- 
consistency, the field equations in brane world models. 
In Sec. Ill, we consider the classical Solar System tests 
in general relativity, namely, the perihelion shift, the light 
deflection and the radar echo delay, for arbitrary spheri- 
cally symmetric spacetimes. In Sec. IV, we analyze the 
classical Solar System tests for the case of various brane 
world vacuum solutions. We conclude and discuss our 
results in Sec. V. 



II. GRAVITATIONAL FIELD EQUATIONS ON 
THE BRANE 



Wc start by considering a 5D spacetime (the bulk), 
with a single 4D brane, on which matter is confined. 
The 4D brane world {'^^^ M , g ^j,,,) is located at a hy- 
persurface [B (X^) = O) in the 5D bulk spacetime 

{^^^M,gAB), where the coordinates are described by 
X"^, A = 0,1, ... ,4. The induced 4D coordinates on the 

brane are .t^, /i = 0, 1, 2, 3. 

The action of the system is given by [2] 



(1) 



where 



Sbulk 



(5)M 



2ki 



d'X, (2) 



and 



Sbr 



(5)c 



(4)M 



+ Lbrane {9a0, V') + -^6 



where fc| 



(3) 

SttGs is the 5D gravitational constant; ^^'R 
and ^^^I/m are the 5D scalar curvature and the matter 
Lagrangian in the bulk, Lbrane {gap, i') is the 4D La- 
grangian, which is given by a generic functional of the 
brane metric ga/B and of the matter fields ip; if ^ is the 
trace of the extrinsic curvature on either side of the brane; 
and A5 and A5 (the constant brane tension) arc the neg- 
ative vacuum energy densities in the bulk and on the 
brane, respectively. 

The energy-momentum tensor of bulk matter fields is 
defined as 



{5)f 



IJ 



IJ 



(4) 



while is the energy-momentum tensor localized on 
the brane and is given by 



(5) 



Thus, the Einstein field equations in the bulk are given 

by [2] 



(6) 



with 

(5) 



Tij = -A^^^^gij + 6{B) [-A6(5)g„ + T„l . (7) 



The delta function 5 {B) denotes the localization of brane 
contribution. In the 5D spacetime a brane is a fixed 
point of the Z2 symmetry. The basic equations on the 
brane are obtained by projections onto the brane world. 
The induced 4D metric is gu = ^^"^gij — ninj, where 
ni is the space-like unit vector field normal to the brane 
hypersurface (*^M. In the following we assume '^^^L^ = 
0. 

Assuming a metric of the form ds"^ = {niUj + 
gij)dx^ dx'^ , with njdx^ = dx the unit normal to the 
X = constant hypcrsurfaces and gu the induced metric 
on X = constant hypersurfaces, the effective 4D gravita- 
tional equation on the brane (the Gauss equation), takes 
the form [2]: 



(8) 



where 5*^^ is the local quadratic energy-momentum cor- 
rection 



Sail — 
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(9) 

and Ej^i, is the non-local effect from the free bulk gravita- 
tional field, the transmitted projection of the bulk Weyl 
tensor Cjajb, Ejj = CjAjsn^n^^ with the property 
Eij Efj^jySjSj as X 0- We have also denoted 
fcl = SttG, with G the usual 4D gravitational constant. 
The 4D cosmological constant. A, and the 4D coupling 
d'^x, constant, k4, are related by A = fc|(A5 + /;;|A^/6)/2 and 
kl = kfXf,/6, respectively. In the limit A^^ — ^ we re- 
cover standard general relativity [2] . 

The Einstein equation in the bulk and the Codazzi 
equation also imply the conservation of the energy- 
momentum tensor of the matter on the brane, Di,T^'^ = 
0, where D,^ denotes the brane covariant derivative. 
Moreover, from the contracted Bianchi identities on the 
brane it follows that the projected Weyl tensor obeys the 
constraint D^Ef," = k^D^^Sf,". 

The symmetry properties of i^^m/ imply that in general 
we can decompose it irreducibly with respect to a chosen 
4- velocity field as 



E, 



u 



(10) 

where k = kr,/k4, h^i, = g^i, -I- u^Uu projects orthogonal 
to w'^, the "dark radiation" term U = —k^E^^u'^u" is 
a scalar, = k^hf^Ea^ is a spatial vector and P^i, = 
—k^ "ft.^,) 1^ — ^hf^^h"^] Eap is a spatial, symmetric 
and trace-free tensor. 
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In the case of the vacuum state we have p = p = 0, 
= 0, and consequently S^i, = 0. Therefore the field 
equation describing a static brane takes the form 



(11) 



where Ci is an arbitrary constant of integration, and we 
denoted 



/ r'^U{r)dr. (19) 

Jq 



with the trace R of the Ricci tensor Rf^i, satisfying the 
condition R = Rif = 4A. 

In the vacuum case i^^^ satisfies the constraint 
DuEfj^ = 0. In an inertial frame at any point on the 
branc we have = 5q and h^^ = diag(0, 1, 1, 1). In a 
static vacuum = Q and the constraint for E^^ takes 
the form [7] 

^D^U+^UA^ + D''P^, + A''P^, = {), (12) 

where A^^ = u"Di,Uij, is the 4-acceleration. In the static 
spherically symmetric case we may choose A^ = A{r)r^ 
and P^j, = P{r) [r^ri, - ^hi^„), where A{r) and P{r) 
(the "dark pressure" although the name dark anisotropic 
stress might be more appropriate) are some scalar func- 
tions of the radial distance r, and is a unit radial 
vector [6]. 

In order to obtain results which arc relevant to the 
Solar System dynamics, in the following we will restrict 
our study to the static and spherically symmetric metric 
given by 



dr^ 



(13) 



where dO^ = dO"^ + s\v?9d4P'. Here 6 and 4> are the 
standard coordinates on the sphere, t G R and r ranges 
over an open interval (rmin, ^max) so that — oo < r^i^ < 
''max < 00. We also assume that the functions ^{r) and 
A(r) are strictly positive and (at least piecewise) differ- 
entiablc on the interval (rmin, rmax)- 

Then the gravitational field equations and the effective 
energy-momentum tensor conservation equation in the 
vacuum take the form [8-10] 



The function Mu is the gravitational mass correspond- 
ing to the dark radiation term (the dark mass). By sub- 
stituting v' given by Eq. (17) into Eq. (15) and with 
the use of Eq. (18) we obtain the following system of dif- 
ferential equations satisfied by the dark radiation term 
U, the dark pressure P and the dark mass Mjj, describ- 
ing the vacuum gravitational field, exterior to a massive 
body, in the brane world model: 

dU _ {2U + P) [2GM + Mu + a{U + 2P) r^] - fAr^ 
dr 



2GM _ Mci „ A^2^ 



^dP 6P 

dr r 



(20) 



and 



dMv 
dr 



= Sar^U, 



(21) 



respectively. 

The system of equations (20) and (21) can be trans- 
formed to an autonomous system of differential equations 
by means of the transformations 

2GM Mu A n o 9. 

r r 6 

p^3ar^P-3r^A, = \n.r. (22) 

We may call /i and p the "reduced" dark radiation and 
pressure, respectively. 

With the use of the new variables given by Eqs. (22), 
Eqs. (20) and (21) become 



— e 



1 



487rG 



U + A, 



(14) 



dq 



(23) 



e-M'^ + 4V4 = ^(/7 + 2P)-A, (15) 



e-Hy"+-^ + - --^—]^^-i^{U-P)-2A, 

(16) 



yl2 yi-y y'X'\ 327rG 



2 ' r 2 J kfXb 



, U' + 2P' 6P 
V = —- 



2U + P r{2U + Py 



(17) 



where we denoted ' = d/dr. 

Equation (14) can immediately be integrated to give 

^_,^^_^_2GM,(.)_A 
r r 3 



and 



d,^_i2, + p)[q+li, + 2p)]_^d._^^^ (24) 
da 1 — q do 

respectively 

Equations (20) and (21), or equivalently, (23) and (24) 
may be called the structure equations of the vacuum on 
the branc. In order to close this system an "equation of 
state" relating the reduced dark radiation and the dark 
pressure terms is needed. Generally, this equation of 
state is given in the form P = P{U). In the new variables 
the functional relation between dark radiation and dark 
pressure takes the formp = 3aexp (29) P [ji/Saexp {26)]. 
We consider several specific solutions below in the con- 
text of the Solar System tests. 
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III. CLASSICAL TESTS OF GENERAL 
RELATIVITY IN ARBITRARY SPHERICALLY 
SYMMETRIC STATIC SPACE-TIMES 

At the level of the Solar System there are three funda- 
mental tests, which can provide important observational 
evidence for general relativity and its generalization, and 
for alternative theories of gravitation in flat space. These 
tests are the perihelion precession of Mercury, the deflec- 
tion of light by the Sun and the radar echo delay observa- 
tions. These tests have been used to successfully test the 
Schwarzschild solution of general relativity. In order to 
constrain the brane world models at the level of the Solar 
System, we first have to study these effects in spherically 
symmetric space-times with arbitrary metrics. In this 
Section, we develop a formalism that can be used for any 
given metric. This useful formalism was first used in the 
Solar System tests applied to a specific vacuum solution 
in Hofava-Lifshitz gravity [27]. 



A. The perihelion precession 

The motion of a test particle in the gravitational field 
on the brane in the metric given by Eq. (13) can be 
derived from the variational principle 

6 J ^ e'^cH'^ - e^r2 _ ^2 (^q2 + gin^ e^^^ds = 0, (25) 

where the dot denotes d/ds. It may be verified that the 
orbit is planar, and hence we can set 9 = tt/2 without any 
loss of generality. Therefore we will use (j> as the angular 
coordinate. Since neither t nor cj) appear explicitly in Eq. 
(25), their conjugate momenta yield constants of motion, 

e'^c^i = E = constant, r^^ = L = constant. (26) 

The constant E is related to energy conservation while 
the constant L is related to angular momentum conser- 
vation. 

The line element Eq. (13) provides the following equa- 
tion of motion for r 



■z I -A 2 12 -\ I V 212 1 

r + e r (p = e (e c t — 1 



(27) 



Substitution of t and ^ from Eqs. (26) yields the fol- 
lowing relationship 



(28) 



The change of variable r = l/u and the substitution 
d/ds = Lu^d/d(p transforms Eq. (28) into the form 



By formally representing e ^ = 1 — f{u), we obtain 

2 



du\' 9 



/("K + ^e-^"'-i2e"'^G(u). 

(30) 

By taking the derivative of the previous equation with 
respect to ^ we find 



where 



F{u) = 



1 dG{u) 

2 du 



(31) 



(32) 



A circular orbit u = uq is given by the root of the 
equation uq = F (uq). Any deviation 5 = u — Uo from a 
circular orbit must satisfy the equation 



d(j)'^ 



dF\ 



5 = O {5^) 



(33) 



which is obtained by substituting u = uq + 5ui Eq. (31). 
Therefore, to first order in 5, the trajectory is given by 



(5 = (5o cos 




(34) 



where and (3 are constants of integration. The angles 
of the perihelia of the orbit arc the angles for which r is 
minimum and hence m or (5 is maximum. Therefore, the 
variation of the orbital angle from one perihelion to the 
next is 



27r 



27r 



(35) 



The quantity a defined by the above equation is called 
the perihelion advance. It represents the rate of advance 
of the perihelion. As the planet advances through (p ra- 
dians in its orbit, its perihelion advances through a4> ra- 
dians. From Eq. (35), a is giv(;n by 



(36) 




or, for small {dF/du}^^^ , by 



"=2 



dF 
du 



(37) 



For a complete rotation we have <p w 27r(l + a), and 
the advance of the perihelion is S(f> = (/) — 2Tr ^ 277(7. In 
order to be able to perform effective calculations of the 
perihelion precession we need to know the expression of L 
as a function of the orbit parameters. Let's consider the 
motion of a planet on a Keplerian ellipse with semi-axis 
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a and b, where b 



and e is the eccentric- 



ity of the orbit. The surface of the elhpse is Trab. 
Since the elementary oriented surface area of the ellipse 
is da = (r X dr) /2, the areolar velocity of the planet is 
\d&/dt\ = |f X dr\ /2 = r2 {d4>/dt) /2 « ^la^^/V^/T , 
where T is the period of the motion, which can be ob- 
tained from Kepler's third law as = An'^a^/GM. In 
the small velocity limit ds w cdt, and the conservation of 
the relativistic angular momentum gives r'^d(l>/dt = cL. 
Therefore we obtain L = 2'kc?\/\ — /cT and = 
(?IGMa (l-e2). 

As a first application of the present formalism we con- 
sider the precession of the perihelion of a planet in the 
Schwarzschild geometry with e" = = 1 — 2GM/ (?r = 
1 - (2GM/c2) u. Hence f(u) = {2GMIc^) u. Since for 
this geometry v + X — we obtain first 



G{u) = 



and then 



2GM 



1 E 



^ , 2GM 



F{u) = 3^u2 + 



GM o GM 

(p-lp- 



(38) 



(39) 



The radius of the circular orbit uq is obtained as the 
solution of the equation 



GM 2 GM 
uo = 3— «o + ^, 

with the only physical solution given by 



(40) 



1± v'l-12G2M2/c4L2 GM , , 
wo = — o ~ (41) 



&GM/c^ 



Therefore 



= TT 



du 



6ttGM 



(42) 



„_„ c^a (1 — e^) ' 

U=Uq V / 

which is the standard general relativistic result. 

B. The deflection of light 

In the absence of external forces a photon follows a 
null geodesic, ds'^ = 0. The afEne parameter along the 
photon's path can be taken as an arbitrary quantity, and 
we denote again by a dot the derivatives with respect to 
the arbitrary afRne parameter. There are two constants 
of motion, the energy E and the angular momentum L, 
given by Eqs. (26). 

The equation of motion of the photon is 



+e" 



= e 



(43) 



which, with the use of the constants of motion can be 
written as 



■2 A-^' 



(44) 



The change of variable r = \/u and the use of the 
conservation equations to eliminate the derivative with 
respect to the affine parameter leads to 



du 
d4> 



\_E^ 



-v-\ 



= P{u) . (45) 



By taking the derivative of the previous equation with 
respect to (j) we find 



where 



(46) 



(47) 



In the lowest approximation, in which the term of the 
right hand side of the equation (46) is neglected, the so- 
lution is a straight line. 



cost; 



(48) 



where R is the distance of the closest approach to the 

mass. In the next approximation Eq. (48) is used on the 
right-hand side of Eq. (46), to give a second order linear 
inhomogeneous equation of the form 



dcjy" 



Q 



cose, 



(49) 



with a general solution given hy u = u {(p) . The light 
ray comes in from infinity at the asymptotic angle = 
— (7r/2 + e) and goes out to infinity at an asymptotic 
angle (p = tt/2 + e. The angle e is obtained as a solution 
of the equation u (7r/2 + e) = 0, and the total deflection 
angle of the light ray is ^ = 2e. 

In the case of the Schwarzschild metric we have v + 
A = and f{u) = {2GM/c^) u, which gives P{u) = 
(2GM/c2) and Q{u) = {3GM/c'^) m^, respectively. In 
the lowest approximation order from Eqs. (48) and (49) 
we obtain the second order linear equation 

d'^u 3GM 2. 3GM 

d^+^=^^os ct>=^{l + cos2ct>), (50) 

with the general solution given by 



cos(/> 3GM /I 



(51) 



By substituting (f) = 7r/2 -|- e, w = into Eq. (51) we 
obtain 



2GM 



(52) 



where we have used the relations cos (7r/2 -|- e) = — sine, 
cos(7r-|-2e) = — cos2£, sine w e and cos2e w 1. The 
total deflection angle in Schwarzschild geometry is 5 = 
2e = AGM/c^R. 
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C. Radar echo delay 

A third Solar System test of general relativity is the 
radar echo delay [28] . The idea of this test is to measure 
the time required for radar signals to travel to an inner 
planet or satellite in two circumstances: a) when the sig- 
nal passes very near the Sun and b) when the ray does 
not go near the Sun. The time of travel of light between 
two planets, situated far away from the Sun, is given by 



J-h 



dy 
c 



(53) 



where h and h are the distances of the planets to the 
Sun. If the light travels close to the Sun, the time travel 

is 



dy, 



(54) 



where v = ce'^" '^^/^ is the speed of light in the presence 
of the gravitational field. The time difference is 

5T = T-To = - r |e[^W-"('-)l/2 - l| dy. (55) 



Since r 
5T=l 



f i?^ , we have 



1 dy. 



(56) 

In the case of the Schwarzschild metric \ = —v, 
exp (A/2 - J//2) = exp {-v) = (l - 2GM/(?rY^ « 1 + 
2GM/c^r, and therefore 



5T ■ 



2GM 



dy 



2GM, JR^ + 11 + I 

— In 



^/W+il-h 

(57) 



Since 



In 



In 



+ q + 



(58) 



where we have used the conditions /l\ ^ 1 and 
R^/q 1, respectively, the time delay is given by 



_ 2GM, 4hh 



(59) 



general spherically symmetric black holes on the brane 
than four-dimensional general relativity. Hence an inter- 
esting consequence of the brane world scenario is in the 
nature of spherically symmetric vacuum solutions to the 
brane gravitational field equations, which could represent 
black holes with properties quite distinct as compared to 
ordinary black holes in four dimensions. Such black holes 
are likely to have very diverse cosmological and astro- 
physical signatures. In the present Section we consider 
the Solar System test properties of several brane world 
black holes, which have been obtained by solving the vac- 
uum gravitational field equations. There are many black 
hole type solutions on the brane, and in the following we 
analyze three particular examples. We will assume, for 
astrophysical constants, the values [28]: 



Mq 
Rq 

c 
G 
a 

e 



1.989 X 10^" kg, 
6.955 X 10® m 
2.998 X 10* m/s. 



6.67 X 10"'' m^kg^^s^^ 
57.91 X 10^ m, 
0.205615 



(60) 



A. The DMPR brane world vacuum solution 

The first brane solution we consider is a solution of the 

vacuum field equations, obtained by Dadhich, Maartens, 
Papadopoulos and Rezania in [6], which represent the 
simplest generalization of the Schwarzschild solution of 
general relativity. We call this type of brane black hole 
as the DMPR black hole. The Solar System tests for the 
DMPR solutions were extensively analyzed in [29], but 
we use the general and novel formalism developed above 
as a consistency check. For this solution the metric tensor 
components are given by 



1 



2m 
r 



(61) 



where Q is the so-called tidal charge parameter. In the 
limit Q — > we recover the usual general relativistic case. 
In terms of the general equations discussed in section II, 
this class of brane world spherical solution is character- 
ized by a equation of state relating dark energy and pres- 
sure: P = —2U. The metric is asymptotically flat, with 
limr^tx) exp [v) = limr_>.oo exp (A) = 1. There are two 
horizons, given by 



IV. SOLAR SYSTEM TESTS OF SPECIFIC 
BRANE WORLD MODELS 

The braneworld description of our universe entails a 
large extra dimension and a fundamental scale of grav- 
ity that might be lower by several orders of magnitude 
compared to the Planck scale [1]. It is known that the 
Einstein field equations in five dimensions admit more 



m 



(62) 



Both horizons lie inside the Schwarzschild horizon = 
2m, < rjj^ < < Ts. In the brane world models there 
is also the possibility of a negative Q < 0, which leads 
to only one horizon r^j^ lying outside the Schwarzschild 
horizon, 



rh+ =m + yjw? + Q > Ts 



(63) 
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In this case the horizon has a greater area than its 
general relativistic counterpart, so that bulk effects act 
to increase the entropy and decrease the temperature, 
and to strengthen the gravitational field outside the black 
hole. 



Perihelion precession 



2. Deflection of light 



We have the following functions 



and 



P{u) = 2mu^ - Qu^ + 



(70) 



For the DMPR black hole we have the following rele- 
vant functions: 



f{u) = 2mu 



(64) 



G{u) = 2mu^ - Qu^ + 



E"^ 1 2mu Qu"" 

- — + — — , (65) 



L2 



i2 



and 



F{u) = 3mu^ - 2Qu^ 



m Qu 



(66) 



respectively. 

uq can be obtained as solution of the algebraic equation 

Stomq - Mo + = 2Qul + -^uo , (67) 

which to first order may be approximated to wq ~ 
GM/{(^L'^), assuming that Q/L"^ < 1. We have used 
the following relationship m = GM/(?. Thus, 5(j) takes 
the form 



5<l) = 



QttGM 



c2a(l-e2) GMa(l-e2) 



2\ ' 



(68) 



which is consistent with the result outlined in [29]. The 
first term in Eq. (68) is the well-known general relativis- 
tic correction term for the perihelion precession, while 
the second term gives the correction due to the non-local 
effects arising from the Weyl curvature in the bulk. 

The observed value of the perihelion precession of the 
planet Mercury is S^pobs = 43.11 ± 0.21 arcsec per cen- 
tury [28]. The general relativistic formula for the pre- 
cession, 5i^GR = 67rGM/c^a (l — e^), gives, using nu- 
merical values given in (60), Si^gr = 42.94 arcsec per 
century. Therefore, the difference /S.ip = Sipots ^SipcR = 
0.17 ±0.21 arcsec per century can be attributed to other 
effects. By assuming that Aip is entirely due to the mod- 
ifications of the general relativistic Schwarzschild geom- 
etry as a result of the five dimensional bulk effects, the 
observational results impose the following general con- 
straint on the bulk tidal parameter Q: 



\Q\< 



GMqu (1 - e^) 



Aip. 



(69) 



With the use of the observational data for Mercury, Eq. 
(69) gives IQI < (5.17±6.39) x 10^ m^, or in the natural 
system of units, with c = h = G = 1, \Q\ < (1.32 ± 
1.63) X 10^0 MeV-^ 



Q{u) = Smu^ - 2Qu^ . 
Solving the differential equation 



(fu 3m o , 2Q o 

_ + «=_cos <^--^cos , 



provides the following exact solution 
cosfli 3m / 1 



(71) 



(72) 
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, „o \ ; cos (h — 7 cos 36 + 36 sin ( 
4i?-^ V4 4 



.(73) 



By substituting (j) = ■jt/2 + e, u = into Eq. (73), 
yields the following total deflection angle 



, , „ 4GM 3nQ 4GM / 3ttQc'^ \ 
5cp = 2e=^-^ = ^(l-^^^j 



c^R 4i?2 c^R 



(74) 

which is in agreement with the results obtained in [30]. 
The best available data on light deflection by the 

Sun come from long baseline radio interferometry [31], 

which gives Sc/^ld = Sc/j^^j^^ (1 + A^d), with A^d < 
0.0002 ± 0.0008, where 6<j)fj^^ = 1.7510 arcsec. There- 
fore light deflection constrains the tidal parameter Q as 
jgj < 16GMRALD/3nc^. By taking for R the value of 
the radius of the Sun in (60), the light deflection gives the 
constraint |Q| < (6.97 ± 27.88) x 10^ m^, or, in natural 
units, 101 < (1.78 ± 7.11) X 10^3 MeV'^. 



3. Radar echo delay 

The delay can be evaluated from the integral in Eq. 
(56), with A = —v, exp (A/2 — 1^/2) = exp (A), so that 



exp (A) 



2GM Q 



2GM Q 



(75) 

so that the time delay, given by Eq. (56), is readily inte- 
grated to yield 



^ 2GM ( ^/w+q + h 



Q_ 
cR 



• (76) 
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Using the approximations /l\ <C 1 and B? /l\ <SC 1, the 
above expression reduces to 



5T'. 



2GM 



In 



i?2 



-1- 



Recently the measurements of the frequency shift of 
radio photons to and from the Cassini spacecraft as they 
passed near the Sun have greatly improved the observa- 
tional constraints on the radio echo delay. For the time 
delay of the signals emitted on Earth, and which graze 
the Sun, one obtains /SiRo = ^t^Ro^ (1 + ^rd), with 
^RD =i (1.1 ± 1.2) X 10-5 [32]. Therefore radar echo de- 
lay constrains the tidal charge of the DMPR brane world 
black hole as \Q\ < 2GMRARDln {Ahh/R'^) /nc^ ~ 
(1.83 ± 1.99) X 10* m2. In natural units we have \Q\ < 
(4.66 ± 5.08) X 10^3 MeV-2. 



B. The CFM solution 

Two families of analytic solutions in the brane world 
model, parameterized by the ADM mass and the 
PPN parameters (3 and 7, and which reduce to the 
Schwarzschild black hole for /? = 1, have been found 
by Casadio, Fabbri and Mazzacurati in [13]. We denote 
the corresponding brane black hole solutions as the CFM 
black holes. 

The first class of solutions is given by 



e'' = l-^, 
r 



(78) 



and 



respectively. The equation of state describing this so- 
lution is t/ = The metric is asymptotically flat. In 
the case 77 > 0, the only singularity in the metric is 
at r = ro = 2m{l + r]), where all the curvature in- 
variants are regular, r = ro is a turning point for all 
physical curves. For 77 < the metric is singular at 
r = Vh = 2m/ (1 — 77) and at ro, with rh > ro, where 
r/i defines the event horizon. We consider this case in 
the analysis outlined below as the CFM2 solution. 



1. Perihelion precession 

The function G{u) in Eq. (30) specializes, for the 
CFMl metric and CFM2 metric in 



Gi{u) 



{l-2mu) (1-|(1-h|7?) mu) 



(1 



-I- 



1 (1 - 2mu) (1 - I (1 



mu) 



i2 



1 



(1 — 2mu) 



(1 - |mu) 



1 



(83) 



and 

G2{U) 



2m(l -I- r])u'^ - 



3 1 — 2m(l -|- r])u 



(1+^r 



i2 

1 - 2m(l -I- ri)u 



{r] + ^Jl- 2m(l + r])u^ 



2 ' 



(84) 



(79) 



respectively, where 77 = 7— l = 2(/3— 1). This solution 
is characterized by an equation of state of the form 



P 

U 



1 



3 m 
4r 



(80) 



3r 



As in the Schwarzschild case the event horizon is located 
at r = r/( = 2m. The solution is asymptotically flat, 
that is lim^^oo = e^°° = lim,.^oo e"^ = e'^°= = 1. We 
consider this case in the analysis outlined below as the 
CFMl solution. 

The second class of solutions corresponding to brane 
world black holes obtained in [13] has the metric tensor 
components given by 



and 



2m ,^ ' 

1 (l+»? 

r 



(81) 



(82) 



respectively. 

The GR limit is obtained for 77 = 0, which gives the ex- 
pected value of the PPN parameter 7=1. The function 
F{u) is quite complicated, but we are only interested in 
the first order corrections to the GR results. Since the 
terms m/L, 77 and E"^ /c^ — 1 are very small, we can ex- 
pand F{u) up to second order, and keep only the term of 
order 0(2) in the product of the infinitesimals. We get, 
for the CFMl and CFM2 solutions, respectively: 



^ ^ ^ m 3(1 + 1) 2 
m 



m ^ 3(1 + 77) ^2 
L'^ m 



Both function can be written as 

7Tll,2 



-^1,2 



-|- 37711, 2M 



(85) 

(86) 

(87) 



where the tilded parameters are defined as: 



mi = 

7Tl2 



(1 



(l-Fr/)m; 



(1 
(1 
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This allows us to use directly the result obtained in GR, 
substituting the values of the tilded parameters. The 
resulting perihelion precession is: 

6-kGM 2nGM 



~ c2a(l - e2) c2a(l - e^) 
6ttGM 6t:GM 



(88) 
(89) 



where we have used the value of the Schwarzschild radius 
of the sun m = GMq/c^, as mentioned above. As in 
the previous section, we consider the case of the planet 
Mercury and assume that the discrepancy between the 
predicted and the observed value is completely due to 
extradimension influence. We thus obtain a limit for the 
value of 77 and consequently for p. In the CFM2 case 
(which is the most constraining) we have the following 
restrictions: 



|77| < 0.004 ±0.005 



(90) 



The T] parameter for the CFMl solution has a bound 
three times larger. 



2. Deflection of light 

The function P{u) in Eq. (45) for the CFMl and the 
CFM2 metric takes the following respective forms: 



_ (1 - 2mu) (1 - I (1 + |?7) mu) 
(1 - fmu) 

[1 - I (1 + |r?) mw] 



P2{u) 



d^L"^ (1 - 
2m(l + r])u^ + 



(91) 



2 

£;2 



1 - 2m{l + ri)u 



{iil+y/l- 2m(l + r?)«) 
The function Q{u) of Eq. (46) can be expanded as: 



(92) 



^ , , -q E'^ m E^ m? 



+3m 



rj 3 E^ 



u^ + 0^{u) (93) 



and 



^ , , E^ m „ E'^m^ 

Q2{u) ^ -»?^;^-3r?^-^M 

77 + 5 777,2 

+37n l + ,-^__ 



7i2+03(y)(94) 



for the CFMl and CFM2 solutions, respectively. 

As stated before, at first order, we can take, for both 
cases, Q{u) = 0, so that the zeroth-order solution is 



cos((/)) 
R 



(95) 



which means the light rays travel on straight lines. As- 
suming, again, that the infinitesimals rj and 777/L are 
comparable, the first order equation is the same obtained 
in standard GR, namely: 



,(1)' 



+ u 



(1) 



8777 



cos (J 



which gives, of course, the standard GR result 



,(1) 



^[3-cos(2</,)]. 



(96) 



(97) 



At second order Eq. (46) can be written as 

,(2)" , „,(2) _ ^2 777 



r] E'^ m 



3 c2 L2 



+ 6 — M^^^TX^^^ +7?- 



,(0)^ 



+ 



37772 



COS 



3 C2 L2 ^3 



cos(0) [2 - cos{2(l))] + m7?^^2^98) 



and 



,(2)' 



,(2) 



E"^ m „ 777 |-n^ i-n 
-^^T^ +6^M(°)lt(l) +377 



777^(0)2 



E"^ m 37772 cos^ 6 , , 

= '1^J2+^ cos(.^) [2 - cos(20)] + 3mr?-^ (99) 

respectively for the CFMl and CFM2 solutions. The 
solutions are provided by 



.(2) 



^ 16i?3 
+21cos(0) +60(/>sin(^) 



16 P 

12cos''(0) 7^— cos2((^) 

3 777 



and 



(2) _ 



Uo = 



rj E"^ 777 
3"^Z2 



R 



2m, 
3.^(100) 



16i?3 

+21cos(0) +60(?isin((/)) 



l2cos'^{4>) - l&rj— cos^U) 

777 



E^ m „ m , , 



respectively. 

The deflection angle is obtained by the equation 

77,(7r/2 + e) = with Scf) = 2e. Assuming L = RE/c 
and substituting again 777 = GMq/c^, we finally find: 



<5</.2 = 4-^ 



, 15 GMq 77 
16 c'R 6 



^ 15 GMq 

^+16'^^ 



(102) 
(103) 



Comparison with data from the long baseline radio in- 
terfcrometry experiments [28], allow us, as in the previ- 
ous section, to place an upper limit for 77 (and thus for 
/3) which is, for the CFM2 solution: 



|?7| < 0.0004 ±0.0016 



(104) 



whilst, as before, the rj parameter for the CFMl solution 
has a bound three times larger. 
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3. Radar echo delay 

The delay can be evaluated from the integral in Eq. 
(56). The integrand, in the case of the CMFl and CFM2 

metrics, is 

exp (^) ~ 1 + (22) + (4 - |r/) m^u^ + 0(^3) , 
exp (^) ~ 1 + (2 + ini)mu + (4 + ^rj) m^u^ + 0{u^) , 

respectively, so that the delay is given by 



and 



„GM, Ahh A T] ^ GM ^ _i 4/iZ 



= 2-^log^ l + - + 2.-^log ^ 



(106) 

These results are to be compared with data from the 
Cassini spacecraft [32] described in the previous section. 
We get a bound for rj in the CFM2 model: 



\r]\ < 0.000021 ± 0.000024 



(107) 



the 77 parameter for the CFMl solution has again a bound 
three times larger. 

So we can conclude that the CFM models add a cor- 
rection to the relevant test quantities which arc directly 
proportional to the PPN parameter that describes the so- 
lution themselves, plus additional terms quadratic in the 
ratio between the Schwarzschild radius and the typical 
dimension of the phenomenon. 



therefore expand F(u) in keeping only terms up to 
second order, which means up to terms 1/c"', and find 



F{u) 



rJ2 

6Gm 



+ ■ 



2 8G2^2 
— u"^. 



(109) 



Therefore uq = F{uo) now becomes a cubic equation for 
uo- The physical solution in up to terms of the order 
{Gm/(?Y is given by 



2Gm 



1 + 



Hence, the perihelion precession is 



dirGm 



1 



4 + 6s - 3s^ 



(110) 



(111) 



0^0^(1 — e^) \ Ss^ 

Comparison with observations yields 

4 + 6s - 3s2 
-0.004 < — ^— ^ < 0.004 , 

3s^ 

which can be achieved by the following parameter choices 



(112) 



-0.528 < s < -0.527, 



or 



2.519 < s < 2.536. 

(113) 

These inequalities place stringent constraints on the pa- 
rameter s, and since s e N, the restrictions strongly vi- 
olate Solar System constrains. Thus, we can safely con- 
clude that the BDM solution is incompatible with Mer- 
cury's perihelion precession. 



2. Deflection of light 



C. The BMD solution 

Several classes of brane world black hole solutions 
have been obtained by Bronnikov, Melnikov and Dehnen 
in [17] (for short the BMD black holes). A particular 
class of these models has the metric given by 



e^= 1- 



2Gm 



2/s 



= 1 - 



2Gm' 



, (108) 



where s G N. The metric is asymptotically flat, and at 
r = rh = IGmjc? these solutions have a double horizon. 

It is also interesting to mention that the DMPR and 
CFM black holes, analyzed above, are part of the fam- 
ilies of solutions classified by Bronnikov, Melnikov and 
Dehnen in [17]. 



1. Perihelion precession 

For the BMD black hole solutions f{u) = {AGm/c^)u- 
{iGm/c^Yv? and F{u) results in an involved expression. 
The equation uo = F{uo) cannot be solved exactly. We 



As above, the functional form of Q{u) is too compli- 
cated to solve the resulting second order differential equa- 
tion. Therefore, we expand Q{u) in 1/c^ and find 



Q{u) 



2Gm 2Gm 6Gm 



/8Ghrf__12G^ 



+ ■ 



8G2m2 
AG^m'' 



uilU) 



where the corresponds to classical general relativity 
and all other terms are contribution from the brane. De- 
spite this complicated form of Q{u) one can solve analyt- 
ically the differential equation 



cose, 



(115) 



which results in a rather lengthy expression. However, by 
evaluating the solution «(</) = 7r/2-t-£) = and by assum- 
ing that e is small, and performing a Taylor expansion, 
the deflection angle becomes 



4GM 



(lA-l)fj 



(116) 
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where the factor up front corresponds to the case of gen- 
eral relativity. Long baseline radio interferomctry con- 
strains the factor in the bracket to be < 1.0017 and there- 
fore s would have to be a very small real number to yield 
compatibility with Solar System tests. This parameter 
choice would also contradict the perihelion precession re- 
sults discussed above. 

Thus, there is no need to further discuss the radar 
echo delay for this solution as it is already incompatible 
with Solar System constraints. Lastly, we would like to 
note that these conclusions also follow by considering the 
PPN parameters for the BMD solution, which are given 
by 7 = 2 and l3 = 2 + 4/5^ - 2/s. Clearly, 7 = 2 violates 
Solar System constrain while the minimum value for /3 is 
given by /3min = 7/4 which is attained when s = 4, again 
incompatible with Solar System tests. 



V. DISCUSSIONS AND FINAL REMARKS 

The study of the classical tests of general relativity pro- 
vides a very powerhil method for constraining the allowed 
parameter space of branc world solutions, and to provide 
a deeper insight into the physical nature and properties 
of the corresponding spacetime metrics. Therefore, this 
opens the possibility of testing braneworld gravity by us- 
ing astronomical and astrophysical observations at the 
Solar System scale. In the present paper, we have de- 
veloped a general formalism that facilitates the analy- 
sis of any given metric and provides the basic theoret- 
ical tools necessary for the in depth comparison of the 
predictions of the brane world models with the observa- 
tional/experimental results. In this context, the classical 
tests of general relativity (perihelion precession, deflec- 
tion of light, and the radar echo delay) were considered 
for specific static and spherically symmetric vacuum so- 
lutions in brane world models. 

First, we analyzed a solution of the vacuum field equa- 
tions, obtained by Dadhich, Maartens, Papadopoulos and 
Rczania (DMPR) [6], which represents the simplest gen- 
eralization of the Schwarzschild solution of general rel- 
ativity. The tightest limit we get on the parameter Q 
came from the perihelion precession of Mercury, and gives 
\Q\ < (1.32 ± 1.63) X lO^o MeV-2, with other bounds 
about three orders of magnitude larger. These results 
represent a significant improvement of the results ob- 
tained in [29], where a method based on the first order 
approximation of the Hamilton- Jacobi result was used. 

Second, we considered two families of analytic solu- 
tions, parameterized by the ADM mass and the PPN 
parameters /3 and 7, found by Casadio, Fabbri and Maz- 
zacurati (CFM) [13]. These solutions can be physically 
seen as "dark pressure-dominated" solutions. Contrary 
to what happens in the DMPR case, here the most strin- 
gent bound came from the radar echo delay data of the 
Cassini experiment, bounding the relevant parameter rj 
to be |r;| < 0.000021 ± 0.000024. 

Thus, these first two models analyzed are not directly 



ruled out by Solar System tests. To get a value that 
is more "physically" meaningful, we can use the match- 
ing conditions between the interior and the exterior of a 
uniform density star in brane world, as derived in [7] to 
obtain the brane tension in terms of the relevant param- 
eter of the different solutions. We get: 



Xb > 



Xh > 



QGMl 
47r|Q|i?| 



27Mqc^ ( 1 



3GMq \ 



87ri?||»7 
9Mqc2 (r?+,/l 



2GMq 



87ri?||eta|(l 



eta) 



(117) 



for the DMPR, the CFMl and the CFM2 solutions 
respectively^ Interestingly enough, the most stringent 
bound comes from the DMPR solution (wc; arc; only inter- 
ested in orders of magnitude here, so we will omit actual 
numerical values and errors): Xb > 10^ GeV^, while the 
CFM solutions gives Ah > lO'* GeV''. Wc do not have a 
definitive explanation for this, but we believe it is related 
to the fact that in the CFM solutions the dark radia- 
tion is basically negligible. Incidentally, the two different 
bounds obtained in the different CFM solutions lead to 
exactly the same bound on the brane tension, thus indi- 
cating that, although causally different, the two solutions 
behave in the same way for what concerns physical phe- 
nomena far away the central source. 

In third place, we considered several classes of brane 
world black hole solutions obtained by Bronnikov, Mel- 
nikov and Dehnen [17] (BMD). The observational Solar 
System test place stringent constraints on the parameter 
space of the model, and these restrictions strongly violate 
Solar System constrains, so we can safely rule out these 
solutions as viable solutions. 

Despite the fact of having a whole plethora of brane 
world vacuum solutions that pass the Solar System ob- 
servational test a few remarks are in order, namely, one 
may basically consider two strategies of obtaining solu- 
tions on the brane. First, the bulk spacetime may be 
given, by solving the full 5-dimensional equations, and 
the geometry of the embedded branc is then deduced. 
Second, due to the complexity of the 5-dimensional equa- 
tions, one may follow the strategy outlined in this pa- 
per, by considering the intrinsic geometry on the brane, 
which encompasses the imprint from the bulk, and con- 
sequently evolve the metric off the brane. In principle, 
the second procedure may provide a well-determined set 
of equations, with the brane setting the boundary data. 
Howcn'er, dc!t(;rmining the bulk geometry proves to be an 
extremely difficult endeavor. 
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